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ABSTRACT 

We examine the variability in the intrinsic absorption in the Seyfert 1 galaxy 
Mrk 279 using three epochs of observations from the Far Ultraviolet Spectro- 
scopic Explorer (FUSE) and two epochs of observations with the Space Telescope 
Imaging Spectrograph on the Hubble Space Telescope. Rather than finding sim- 
ple photoionization responses of the absorbing gas to changes in the underlying 
continuum, the observed changes in the absorption profiles can be understood 
more clearly if the effective covering fraction of the gas in all emission compo- 
nents, continuum and broad and intermediate velocity width emission lines, is 
accounted for. While we do not uniquely solve for all of these separate covering 
fractions and the ionic column densities using the spectral data, we examine the 
parameter space using previously well-constrained solutions for continuum and 
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single emission component covering fractions. Assuming full coverage of the con- 
tinuum, we find that of the two velocity components of the Mrk 279 absorption 
most likely associated with its outflow, one likely has zero coverage of the inter- 
mediate line region while the other does not. For each component, however, the 
broad line region is more fully covered than the intermediate line region. Changes 
in the O VI column densities are unconstrained due to saturation, but we show 
that small changes in the nonsaturated C IV and N V column densities are con- 
sistent with the outflow gas having zero or partial covering of the intermediate 
line region and an ionization parameter changing from ~0.01 to ~0.1 from 2002 
to 2003 as the UV continuum flux increased by a factor of ~8. The absence of a 
change in the C III absorbing column density is attributed to this species arising 
outside the Mrk279 outflow. 

Subject headings: galaxies: active — galaxies: individual (Mrk 279) — galaxies: 
Seyfert — quasars: absorption lines — ultraviolet: galaxies — X-ray: galaxies 

1. Introduction 

Mass outflows from low-rcdshift active galactic nuclei (AGNs) are inferred from blueshifted, 
variable X-ray and UV absorption which often shows evidence for partial coverage of the 
contintuum source in nearly one half of all nearby Seyfert galaxies (see review in Crenshaw, 
Kraemer, & George 2003 and recent results in Dunn et al. 2007). In addition to possibly 
playing a strong role in the accretion process (Blandford & Begelman 1999, 2004), AGN 
outflows likely have substantial impacts on their environments, providing feedback energy to 
regulate the formation of galaxies and clusters of galaxies (Granato et al. 2004, Scannapieco 
& Oh 2004, Scannapieco, Silk, & Bouwcns 2007); to set the luminosity-temperature relation 
in galaxy clusters (Cavaliere et al. 2002); to distribute metals in galaxy clusters (Moll et al. 
2007) and into the intergalactic medium (Adelberger et al. 2003). 

Theoretical models of AGN outflows, from accretion disk winds (Konigl & Kartje 1994, 
Murray et al. 1995, Proga 2000, 2003, Proga, Stone, & Kallman 2000, Proga & Kallman 
2004), or from ablation from the obscuring torus (Krolik & Kriss 1995, 2001) place the ab- 
sorbing material at different distances from the central engine and with different orientations 
with respect to the broad line region (BLR). A large sample of high resolution X-ray and 
UV spectra of onflows in nearby Seyferts from Chandra, XMM-Newton, the Far Ultraviolet 
Spectroscopic Explorer (FUSE) and from the Space Telescope Imaging Spectrograph (STIS) 
onboard the Hubble Space Telescope has been building in the literature over the last decade, 
providing much insight into the physical conditions in these onflows. However, a broad un- 
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derstanding of these outflows has yet to emerge, and several outstanding questions remain, 
namely: What is the relationship between the X-ray- and UV- absorbing gas?; What is 
the location and structure of the absorbing material?; Is the observed variability in some 
absorption profiles a response to flux variations in the central engine (Krolik & Kriss 1995, 
Shields & Hamann 1997, Crenshaw et al. 2000, Kraemer et al. 2002) or to bulk motions of 
the absorbing gas (Crenshaw & Kraemer 1999, Kraemer et al. 2001, Gabel et al. 2003b)? 

In general, it is difficult to associate a particular UV absorption component with the X- 
ray absorber unambiguously, in part due to relatively low spectral resolution in X-ray data 
and in part because the ionization modeling is sensitive to the shape of the UV-to- X-ray 
continuum (Kaspi et al. 2001). However, with the exception of campaigns on Mrk 279, using 
FUSE, STIS, and Chandra vnih. both the High Energy Transmission Grating (HETG) (Scott 
et al. 2004, Gabel et al. 2005a, S04 and G05 hereafter) and the Low Energy Transmission 
Grating (Costantini et al. 2007), on NGC 7469 (Scott et al. 2005) using FUSE, STIS, and 
Chandra/ HETG. the long-term observing campaign on NGC 3783, using FUSE, STIS, and 
Chandra/ HETG (Kaspi et al. 2002, Netzer et al. 2003, Gabel et al. 2003a, 2005b), and 
the campaign on NGC 4151 using FUSE, STIS, and Chandra/ HETG (Kraemer et al. 2005, 
Kraemer et al. 2006, Crenshaw & Kraemer 2007), UV and X-ray observations are generally 
obtained at different times. 

Several observations of the absorption proflles taken over a period of time are particularly 
useful for addressing the last two questions mentioned above, as they provide contraints on 
the response time of the gas to variability in the central engine and thus limits on the distance 
of the absorber from the continuum-emitting region. Mrk 279 {zem = 0.0305) has been the 
subject of an intensive campaign of simultaneous observations in the UV and X-ray (S04, 
Kaastra et al. 2004, Arav et al. 2005, 2007, G05, Costantini et al. 2007). In this paper, we 
use several epochs of UV spectral obsevations of Mrk 279 over 1999 to 2003 from FUSE and 
STIS presented in S04 and G05 to investigate the variabihty of the column density and/or 
covering fraction of the outflow over time. We choose to concentrate on the UV data because 
the 2002 epoch of Chandra observations coincided with a low flux state. The grating data 
therefore have signal-to-noise (S/N) too low to perform detailed analyses (S04). 

The spectral coverage of FUSE provides information on H I absorption in Ly f3 and 
higher order Lyman transitions, and C III A977, while the STIS observations cover Lya 
and the N V, Si IV, and C IV doublets. The different epochs permit us to investigate the 
effects of UV continuum variability on the intrinsic absorption. The analysis presented here 
will also take into account changes in absorption profiles that may be a result of differing 
contributions of the various emission sources, i.e. continuum versus broad lines, in different 
observation epochs. 
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2. Continuum and Emission Lines 

Tabled] summarizes the data used in this analysis. See S04 and G05 for full presentations 
of the FUSE and STIS spectra. The 1999 FUSE spectrum and the 2002 and 2003 FUSE + 
STIS spectra are shown in Figure [H 

We fit a power law continuum of the form fx = /looo ( — to the 2003 FUSE and 

V 1000 A / 

STIS spectra simultaneously, where the normalization, /looo; is the flux at 1000 A, and we 
include the Galactic extinction law of Cardelli, Clayton & Mathis (1989) with Ry = 3.1 and 
E ( B - V ) = 0.016 (Schlegel, Finkbeiner, & Davis 1998) in the continuum fits. The best fit 
continuum slope and normalization are a = 1.261 ± 0.005 and /looo = 1-230 ± 0.002 x 10~^^ 
ergs s~^ cm~^ A~^. For the 1999- and 2002-epoch data we found (a, /looo) = (1-60 ± 
0.02, 1.325 ± 0.004 x lO^^^) and (0.86 ± 0.02,0.133 ± 0.001 x lO^^^), respectively. These 
results are shown graphically in Figure O The 1999 and 2003 epochs are high flux states, 
and the continuum becomes bluer as it brightens. 

We fit Gaussian profiles to the emission lines in the combined FUSE+STIS spectrum. 
The emission lines we identify and fit are listed in Table [2l See Table 3 of S04 for the fits to the 
1999- and 2002-epoch spectra. An F-test confirms that the most prominent features require 
contributions from both a broad velocity component (BLR, FWHM~ 7500 — 10000 km s^^) 
and an intermediate one (ILR, FWHM~ 3000 km s"^). The resulting values of the reduced 
chi-squared, xl = 1-6 and 1.5, for 2002 and 2003 data, respectively, are significantly improved 
from the values of ~2 achieved with no ILR component. We do not attempt to add any more 
components to the fit to improve the xl because there are no remaining systematics in the 
fit residuals and because these added components would not be well-motivated physically. 
The full spectra with fits are shown in Figure [H 

We note that some of the fitted emission lines show velocity shifts of hundreds of km s^^ 
from their 2002 positions. We performed an emission line fit for the 2003 data with the line 
centroids tied to their 2002 values. This led to a larger xl in the final fit (1.6 versus 1.5), 
but did not change the results of the absorption line analysis presented below. Since we will 
treat our emission model as phenomenological and because S04 did not constrain the line 
centroids in this way (e.g. in fitting the 1999, 2000, and 2002 epoch spectra), we will use 
the emission model with the emission line velocity shifts in place. 

We show the emission profiles for the 2002 and 2003 spectra in Figures [3ll6l For all 
these ions, the continuum flux is clearly stronger relative to the emission lines in the high 
flux states in 1999 and 2003. The O VI intermediate line emission is stronger relative to the 
broad line in 2002 than in 1999 or in 2003, when the continuum fluxes were higher. In fact, 
the ILR component is very weak in the 2003 spectrum. The same is true for the Lya and 
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N V intermediate lines. They too are stronger relative to their respective broad lines in the 
2002 low flux state than in 2003, the only two epochs we can compare with STIS data. For 
C IV, the difference between the low and high flux states is less pronounced. Figure El shows 
that the intermediate width line is only slightly more prominent relative to the broad line 
in 2002 than in 2003. 



S04 identified several velocity components to the absorption in the Mrk 279 spectrum 
extending over -600 km s~^ - +100 km s~^ with respect to the systemic redshift of the AGN, 
taken to be 2 = 0.0305 ± 0.0003. These components, based on the velocity structure in the 
Lj/3 profile, were labeled 1, 2, 2a, 2b, 3, 4, 4a, and 5. G05 added a component 2c to this, 
based on additional velocity structure in the C IV and N V absorption. S04 concluded that 
Components 1 and 4 {v ^ +90 km s~^ and —450 km s~^ in Ly/?, respectively) are likely 
associated with the disk and possibly the halo of the host galaxy of the AGN and/or its 
companion, MCG +12-13-024, while the origins of Components 3 and 5 (v ~ —385 km s~^ 
and —540 km s~^ in Ljf3, respectively) are ambiguous. The absorption from highly ionized 
species over the entire velocity range is taken to be intrinsic by G05, as is the H I absorption 
that is uncontaminated by any possible contribution from low-ionization components seen 
in C II, C III, N III, Si II, and Si III, i.e. v ^ 200 - 300 km s'^. We will follow that 
convention in this paper, referring to Components 2 and 4 in the highly ionized species to 
be the absorption complexes centered at v ~ —300 and v ~ —480 respectively. 

The most general form of the intrinsic absorption profile for a particular species accounts 
not only for partial covering, but also for multiple emission sources: 



where / is the line intensity at velocity bin j, R is the ratio of source i to the total flux at 
j, r is the optical depth of the absorber at j, and C is the effective covering fraction of the 
absorber with respect to source i at velocity bin j (S04, G05). 

The effective covering fraction at bin j is the weighted sum of individual covering frac- 
tions summed over all emission sources (G05): 



G05 found that changes in the Lya absorption profile between the 2002 and 2003 epochs 
could be explained by the change in the effective covering fraction arising from the change 



3. Intrinsic Absorption 



/, = 5^i?}(qe-^ + i-q) 



(1) 




(2) 
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in the relative contribution of the continuum and broad hne region to the total flux in each 
epoch. No changes in the individual emission component covering fractions or in the H I 
column density are necessary to produce the observed reduction in depth of the absorption 
trough in Components 2-2a between the 2002 and 2003 epochs. 

However, we found that this method cannot account for the changes in the absorption 
profiles of O VI, N V, and C IV between the 2002 and 2003 epochs if only the continuum and 
emission line sources are considered separately. Below we will examine this ionic absorption 
variability while treating the intermediate width emission line components as a separate 
emission source. For O VI, we also include the 1999 epoch for which we have additional 
FUSE data. We will use the well-constrained covering fraction and column density solutions 
derived by G05 from the high flux state 2003 data shown in Figure [71 These solutions are 
the global CNO flts, the simultaneous flts to the O VI, N V, and C IV doublets described 
by G05. The continuum covering fraction is effectively unity over the entire velocity range 
of the absorption while the total emission line covering fraction varies. For the remainder of 
this paper, however, we shall set C'^°'^^ = 1 (see also Arav et al. 2007). 

First, however, in the following section we will describe the general observed variation of 
the intrinsic absorption. For each species, in each epoch, we show and discuss the normalized 
absorption proflles in the following flgures and sections under the assumption C^^^ = C^^^ = 
(-Ycont _ know this assumption is not strictly true, given the covering fraction and 

column density solutions shown in Figure [3 but we use this as a heuristic model and a 
starting point for further analysis. In Section HJ we discuss photoionization models for this 
outflow, and in Section [5] we discuss the importance of these photoionization effects relative 
to changes in the effective covering fraction of the absorption resulting from variations in the 
underlying emission sources, primarily the broad and intermediate velocity width emission 
lines. 



3.1. Description of Normalized Absorption Profiles 

3.1.1. VI 

Figure [8] shows the O VI proflle. Given the equal depths of the red and blue doublet 
components, it is clear that the O VI absorption in the core of Component 4 is completely 
saturated in the 2002 and 2003 epochs, and less so in the 1999 high flux state. If C^^^ = 1 
as assumed here, the increase in absorption depth of the saturated core between 2002 and 
2003 indicates the total effective covering fraction increased. This is a contradiction to our 
assumption here, i.e. unity covering fractions for all emission components for all epochs. 
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confirming what we know from Figure U\ namely that C^^^ and possibly also C^^^ are not 
unity. 

The O VI absorption is saturated in Component 2 in all three epochs. For C^^^ = 1, 
the top panel of Figure [8] shows that the absorption depth is largest in 2003 and smallest in 
2002. This is true of Component 4 as well. These are the epochs of the smallest and largest 
ILR/BLR flux, respectively, indicating that the profile change is caused by the change in 
this emission line ratio, contributing to a change in the total effective covering fraction. 

3.1.2. NV 

In Component 4, the entire N V profile is marginally saturated in 2002 and unsaturated 
in 2003. For C™°* = C^^r ^ ^-^ilr = 1^5 shown in Figure [H the absorption depth is 
similar but nonzero in both epochs. This indicates once again that the emission line covering 
fractions are nonunity, but also that the change in the effective line covering fraction between 
the two epochs is not as dramatic as for O VI. 

For Component 2, the N V profile is saturated in 2002, but in 2003 it is saturated in the 
core and in the red wing of the profile, while the blue wing is unsaturated. The top panels 
of Figure [9] show the absorption depth in 2003 is greater than in 2002, particularly in the 
consistently saturated regions. This could be the result of the changing ILR/BLR flux ratio 
and the resulting variation in the effective line covering fraction, as described above for the 
O VI profile. 

3.1.3. CIV 

In 2002, Component 4 of the C IV profile is possibly saturated in the red wing and un- 
saturated elsewhere, while in 2003 it is unsaturated over most of the profile. The absoprtion 
depth decreases from 2002 to 2003, indicating that the C IV column density decreased as 
the continuum flux increased by a factor of ~ 8 between 2002 and 2003. The profile depth 
change between 2002 and 2003 is less than expected for a simple column density response 
to this fiux increase, leading to the conclusion that the total effective covering fraction de- 
creased over this time. Although there is only a small change in the ILR/BLR emission flux 
ratio (Fig. [6]) the continuum level relative to the emission lines rises signiflcantly between 
2002 and 2003. The effective covering fraction thus rises since C*^""^* = 1. 

In Component 2, the C IV proflle is unsaturated in both 2002 and 2003, and the ab- 
sorption depth is unchanged between these two epochs, most likely also due to ionization 
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changes balancing an increase in the effective covering fraction. 

3.1.4- Lyman series 

For H I, the only part of the absorption complex that can be reliably attributed to 
intrinsic gas is Component 2 {v ^ 200 — 300 km s~^). In this velocity component, the Lya 
absorption depth is noticeably weaker in 2002 than in 2003. As pointed out by G05, this is 
explained by a lower effective covering fraction in 2002 due to a larger relative contribution 
by the BLR relative to the continuum and the fact that (7™°* = 1 while C^^^ < 1. For our 
assumption C^^^ = C^^^ = 1 ^ shown in Figure [TT], we see this reported change in the Lya 
profile, and little change in the Ly/3 profile, which we would expect to be more influenced 
by the changes in the O VI emission components. 

For other lines further along the Lyman series from Lj(3, there is little or no underlying 
emission line flux beneath these profiles. Examining Ly7, a Lyman series line unlikely to be 
affected by saturation, we find no discernable difference in the profiles from 1999 to 2002, 
but the profile does appear less deep in 2003. However there is no difference in LjS or Lye 
between the two high flux states, 1999 and 2003. 

4. Photoionization 

To investigate a wider range of possiblities for the various emission source covering frac- 
tions, and to disentangle the effects of photoionization from those of covering fraction and 
emission source variation on the various species considered above, we computed photoionza- 
tion models similar to those of Krolik & Kriss (1995, 2001) using the XSTAR code (Kallman 
2000). We calculate the photoionization state of the gas using the supersolar abundance 
found by Arav et al. (2007) spectral energy distributions (SEDs) appropriate to the 2002 
and 2003 FUSE and Chandra observations shown in Figure [T2l Both SEDs are normalized to 
the measured FUSE and Chandra spectra over the spectral range of these instruments. The 
parameters of the F spectra are discussed in Section [21 We connect the 2002 FUSE and 
Chandra measurements with auv-x-ray = 0.57; and we approximate the 2002 Chandra dou- 
ble power law spectrum (S04) with a single power law with a = 1.47 from 0.5 to 50 keV and 
we extend the SED to ~150 keV with a = 0. We approximate the 2003 Chandra power law 
+ modified blackbody continuum (Costantini et al. 2007) with a smooth power law: a = 0.5 
connecting the FUSE and Chandra data to 0.7 keV and a = 0.88 to ~150 keV (Costantini 
et al. 2007, Arav et al. 2007). For energies lower than the STIS bandpass and higher than 
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the Chandra bandpass, we use cut-offs for a standard AGN spectrum (cf. Costantini et al. 
2007). We executed tfie model witli grids in input column density and ionization parameter 
from logA^ = 18 to logA^ = 21 and from U = 0.01 to U = 10, respectively. 

In Figure [131 ^6 show results from the log = 20 run. We plot the ratio of the column 
densities of C III, C IV, N V, and O VI in the XSTAR simulations using SED2002 to the 
values from the XSTAR simulations using SE1I)2003 

function of the ionization parameter 
in the SED2002 simulations. To take this ratio, we compare the column densities in each 
simulation at equivalent ionization parameters given the different SEDs. The number of 
ionizing photons supplied by SED2002 is about one-tenth the number supplied by SED2oo3- 
Therefore, we take the ratio of the column densities of each ion at [/2002 and [/2003 = 10?72oo2- 
These are plotted against t/2002 in Figure O 



5. Effective Covering Fraction versus Photoionization 

We now seek to establish whether the changes in the absorption profiles of the highly 
ionized species C IV, N V, and O VI between the 2002 and 2003 epochs can be explained by 
partial covering of some portion of the line emission combined with changes in the relative 
contributions of the component broad and intermediate line fluxes to the total emission 
profile, or if some change in the column density in response to changes in the continuum flux 
is required by the data. To address this, we use the well-constrained 2003 covering fractions 
and column densities derived by G05 and attempt to reproduce the 2002 absorption profiles 
in two ways: (1) by allowing the covering fraction of the ILR to be different from that of the 
BLR; and (2) by varying the column densities according to changes in the ionizing flux as 
predicted by our photoionization models. From Figure [T31 we see how we expect the column 
densities of each ion to change under photoionization alone. 

From the emission line fits, we consider three emission sources, the continuum, the BLR, 
and the ILR. Thus, Equation [T] becomes 

We assume (7'^°'^* = 1, consistent with G05. For doublets, however, the resulting system of 
equations is still underconstrained. Therefore, we assume three possibilities for C^^^: (1) 

^ILR _ ^BLR. C^LR = O.SC^LR. g^^d (3) C^^R = Q. 

The overall goal of this experiment is to reproduce the 2002 spectrum from the column 
densities derived from the higher quality 2003 data. We attempt to do this by varying both 
the changes in the ionic column densities using our photoionization modeling as a guide, and 
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by varying the covering fractions of the intermediate hne region using the three cases above. 
Here is a summary of our method: 

• We begin with the well-constrained solutions for the the emission line covering fractions 
from the 2003 data presented by G05. In that solution, the line covering fractions of 
C IV, N V, and O VI were tied to the same value, 0.64 and 0.76 for the centers of 
Components 2 and 4, respectively. 

• We use the new emission line fits to the 2003 data presented in Section O including 
intermediate line components, to derive C^^^ and C^^^ from the values of R^^^ and 
R^^^ at each velocity bin, from the total line covering fraction given by G05. Note that 
in doing this, we are formally calculating different covering fractions for C IV, N V, 
and O VI because the emission line flux underlying each set of doublets is different. 
We will return to this point in the discussion in Section [Hi 

• We apply these to the 2002 epoch given the values of R^^^ and R^^^ derived from that 
spectrum to examine how the profiles may have changed from changes in these ratios 
alone. 

• We also examine how the profiles change if the column densities of the various ions 
have changed as well, by reasonable factors gleaned from our photoionization modeling 
as shown in Figure [T31 

• We center the x-axis of Figure [T3l and concentrate our attention on the case log(t/2003) ~ 
— 1 (and Iog(f72oo2) ~ —2) since this is the result of the detailed photoionization mod- 
eling and abundance analysis of Arav et al. (2007). 

6. Discussion 

The results for O VI, N V, and C IV are shown in Figures [T^lfT6l In each lower left panel, 
we show the blue and red components of the doublet profiles from the 2002 data. In the other 
three quadrants, we show the blue component of the 2002 profile in black. Superimposed on 
this are the attempts to reproduce this profile from the parameter constraints from the high 
signal-to-noise 2003 data described above, under variations in both the column densities and 
the covering fractions. The column density changes from 2002 to 2003 increase clockwise 
from top left to bottom right, with one panel for each ion showing how the profiles would 
change if no change in the ion column density occured between 2002 and 2003. 

Figure [H] shows that because the O VI profile is saturated, the column density is un- 
constrained. For Component 2, regardless of the column density the only value of C^^^ 
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consistent with the profile is zero, while for Component 4, C^^^ = 0.5C^^^. The photoion- 
ization model predicts an increase of about a factor of four in the O VI column density from 
2002 to 2003, given the change in the continuum level over that period. 

For N V, shown in Figure [151 the best fit solution for both Components 2 and 4 is a 
column density decrease of ~30% from 2002 to 2003 and C^^^ = 0, although C^^'^ = O.SC^lr 
also provides an acceptable fit for Component 4. Our pure photoionization model would 
predict a decrease of about one half that amount in the N V column density between 2002 
and 2003. 

For C IV, Figure [16] shows that a decrease in the column density from 2002 to 2003 is 
required to match the profiles in both components, regardless of the value of C^^^. A change 
in the column density by larger than a factor of ~2, however, is not consistent with the 
observed absorption depths. Our photoionization modeling indicates that we would expect a 
larger decrease in the C IV column density from photoionization alone, a factor of nearly six. 
This covering fraction is relatively poorly constrained because the relative contributions of 
the BLR and ILR change very little between the 2002 and 2003 epochs, as noted above. Of 
our models, the best fits to the 2002 spectrum come from A^(C IV)2oo2=2 x A^(C IV)2oo3 and 
C^LR _ g Component 2 and C^^^ = C^^^ for Component 4. This ILR covering fraction 
solution for Component 4, however, is not significantly different from C^^^ = 0.5C^^^. 

In the formal solution for these covering fractions we assumed that the effective line 
covering fractions for O VI, N V, and C IV are all the same in the 2003 epoch but that the 
resulting individual BLR and ILR covering fractions are not necessarily equal, and that those 
individual BLR and ILR covering fractions did not change between 2002 and 2003. However, 
from the success of the G05 global CNO covering fraction fits, the covering fractions for each 
of these highly ionized species are likely to be comparable. Thus, we summarize the results 
detailed above in this context and find that (1) the BLR covering fraction of C IV, N V, 
and O VI in the cores of Components 2 and 4 are ~0.8-0.9; (2) the ILR covering fractions 
are zero for Component 2 and ~0.4 for Component 4; and (3) the column density in O VI 
decreased by an indetermiate factor while the column density in N V decreased by one third, 
and that of C IV decreased by a factor of ~2 between the two epochs. Should we choose 
to constrain the individual broad and intermediate emission line covering fractions of C IV, 
N V, and O VI to be equal (like G05, but unlike S04), we could use this set of three doublets 
to solve for six variables, the continuum, broad line, and emission line covering fractions and 
the column densities of O VI, N V, and C IV as a function of velocity across the absorption 
profiles. We leave this for future work. 

The C III profile presents somewhat of a puzzle in light of the results we have outlined 
thusfar. We expect the column density of this species to be quite sensitive to the factor of 
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~8 increase in the ionizing flux between the 2002 and 2003 observations. Figure [13] shows 
that for f/2002 = 0.01, the C III column density should decrease by a factor of about 25-30 
from 2002 to 2003. In Figure [T71 we demonstrate that we see no significant change in the 
C III profile in either Component 4 or in Component 2. 

One way to resolve this problem is to invoke a long photoionization response time for 
the gas. To test this possibility, we follow an equilibrium photoionization model of the gas 
over time. We began with a model that matches the well-constrained 2003 column densities 
and dropped the photoionizing flux by a factor of 8 to evolve the absorber back in time to the 
2002 epoch. Following the trend of ionic abundances, we find that the C III column density 
should decrease by a factor of two in about one hour for Tig = 3 x 10"^ cm~'^, in about one 
day for = 1500 cm~^ and in about one year for = 5 cm~^. Therefore, for a very small 
change in the C III column density, a small total density is required if the absorbing gas is 
in equilibrium. This density, along with logU = — 1 and the luminosity of Mrk 279, require 
that the gas lie at ~2 kpc from the ionizing source if it is part of the Mrk 279 outflow. This 
would certainly present us with a conundrum given the nonunity covering fractions we have 
consistently found for the emission lines in absorption Components 2 and 4. 

The stable C III profiles are not likely solely the result of problems with the normaliza- 
tion of the spectra. We confirm this by comparing the profiles of two nearby ISM lines in 
the three FUSE epochs. In Figure [18] we show two prominent but unsaturated ISM lines, 
S III A1012 and Ar I A1048, in all three FUSE observation epochs. We do see some normal- 
ization problems, particularly in the Ar I region, likely attributable mainly to the low S/N 
in the 2002 data. The absorption depth in the 2002 Ar I profile is higher than in the other 
epochs when it is not expected to vary at all. This may be a contributor to the apparent 
lack of variability in the C III profiles. We note, however, that the deepest absorption in the 
C III profile peaks at -540 km s~^, -460 km s~^, and at -385 km s^^. As discussed by S04, 
the low ionization absorption at these velocities may not be part of the Mrk 279 outflow but 
may in fact be attributable to either the host galaxy of the AGN or to a companion galaxy, 
MCG +12-13-024. 

Photoionization models of quasar line emission indicate that intermediate line region 
lies at approximately 1 pc from the central ionizing source, a factor of ~10 further than than 
the broad line region (Brotherton et al. 1994). However, the UV and X-ray observations of 
the Seyfert galaxy NGC 4151 place the D+E subcomponents of its outflow absorption in 
an accretion disk wind which is the ILR and which lies at only 0.1 pc from the continuum 
source (Kraemer et al. 2006, Crenshaw & Kraemer 2007, see also Corbin 1997). For Mrk 279, 
we have found ILR covering fractions that are consistently smaller than those of the BLR 
and continuum, and in some cases C^^^ = 0. This would place the absorbing gas within or 
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interior to the ILR. Assuming that the width of the intermediate hnes, ~2900 km s~^, is 
produced gravitationally, a black hole mass of 10^'^^ Mq (Peterson et al. 2004) gives a radius 
of ~0.02 pc for the ILR, an estimate that is only slightly larger than measured constraints on 
the size of the BLR, ~0.01 pc (Maoz et al. 1990). We note, however, that this estimate of the 
ILR size is influenced by the relatively broad ILR emission features, certainly as compared 
with NGC 4151. 



7. Summary 

We find that some of the variability, or lack thereof, in the absorption profiles of the 
various species present in the FUSE and STIS data for Mrk 279 can be understood in terms 
of changes in the effective covering fraction when all emission components are included in 
the absorption model. Accounting for individual covering fractions in the continuum, broad 
lines, and intermediate velocity width lines, we find less dramatic changes in the column 
densities of the highly ionized species are required by the absorption data than those that are 
predicted by photoionization modeling alone. In particular, the contribution of the covering 
fraction of the intermediate velocity width emission lines to the overall effective covering 
fraction for a given ion and the changes in the relative contribution of this component to the 
overall emission between the two epochs can account for a significant portion of the changes 
observed in the absorption profiles. 

If the C III absorption is intrinsic to the outflow, the UV continuum flux variation 
between 2002 and 2003 is not commensurate with an absorber of very high ionization pa- 
rameter or of very large density, placing this absorbing structure at a substantial distance 
from the ionizing source. We consider this unlikely, however, given the nonunity broad and 
intermediate emission line covering fractions derived for the intrinsic absorption. 

J. E. S. acknowledges the support of a National Research Council Associateship held 
at NASA Goddard Space Flight Center and the support of the Jess and Mildred Fisher 
Endowed Chair in the Biological and Physical Sciences, held at Towson University. 
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Table 1. Observations of MRK 279 



Instrument 


ID 


Start Date 


Exp. (s) 


FUSE 


P1080303 


1999-12-28 


61139 


FUSE 


C0900201 


2002-05-18 


47414 


FUSE 


D1540101 


2003-05-12 


91040 


STIS/E140M 


O6JM01 


2002-05-18 


13193 


STIS/E140M 


O8K10 


2003-05-13 


41386 
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Table 2. Emission Line Fits to 2003 FUSE and STIS Spectra of 

Mrk 279 



Line 


\ 


r lux 


velocity 


r VV rliVl 








(km s ^) 


(km s ^) 


S VI 


933.378 


4.09±1.64 


683±1782 


7569±192 


S VI 


944.523 


2.04±0.82 


683±1782 


7569±192 


C III 


977.020 


33.55±2.82 


384±458 


9206±146 


N III 


989.799 


8.06±2.52 


417±1066 


9206±146 


Ly/3 broad 


1025.7223 


31.52±2.47 


-79±7 


8731±61 


Ly/3 int. 


1025.7223 


5.77±0.67 


-469±11 


3208±22 


Ly/? narrow^ 


1025.7223 


0.56±0.33 


145±22 


697 


VI broad 


1031.9265 


71.66±2.58 


475±25 


7569±192 


VI int. 


1031.9265 


10.29±1.25 


-193±100 


2783±97 


VI broad 


1037.6155 


35.83±1.29 


475±25 


7569±192 


VI int. 


1037.6155 


5.14±0.62 


-193±100 


2783±97 


S IV int. 


1062.66 


11.12±0.31 


-416±132 


3274±560 


S IV int. 


1072.97 


11.12±0.31 


-416±132 


3274±560 


He II 


1085.15 


28.52±0.79 


-213±101 


7569±192 


Si II 


1192.3 


3.57±0.79 


446±200 


9385±269 


Lya broad 


1215.6701 


307.69±2.38 


-79±7 


8731±61 


Lyo; int. 


1215.6701 


120.23±0.52 


-469±11 


3208±22 


Lya narrrow^ 


1215.6701 


7.09±0.48 


145±22 


697 


N V broad 


1240.40 


106.66±2.80 


-56±53 


9908±19 


N V int. 


1240.40 


5.84±0.89 


-15±215 


2956±72 


Si II 


1260.42 


5.32±0.61 


446±200 


9385±269 


Si II+O I 


1304.35 


9.82±0.89 


-69±140 


4020±732 


CII 


1335.30 


1.68±0.59 


-219±352 


2422±207 


Si IV 


1393.76 


24.77±0.30 


-424±6 


9908±19 


Si IV 


1402.77 


12.38±0.15 


-424±6 


9908±19 


OIV] 


1402.06 


21.93±0.17 


-600±67 


3774±47 


C IV broad 


1549.05 


30.21±0.79 


4±19 


9908±19 


C IV int. 


1549.05 


8.87±1.21 


-304±12 


2956±72 


He II 


1640.5 


80.96±3.30 


-11±122 


9909±19 
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^Flux in units of 10 ergs cm ^ s ^ 

^Velocity relative to systemic redshift, z — 0.0305 

^FWHM fixed to 2002 May value 
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Fig. 1.— 1999 FUSE spectrum (top) and 2002 and 2003 FUSE + STIS spectra (middle and 
bottom) of Mrk 279 with errors (dotted line) and continuum and emission line fits (red) 
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Fig. 2. — Variation in continuum flux and UV spectral index 1999-2003 
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Fig. 3. — O VI emission line profile in 1999, 2002, and 2003 FUSE spectra (cyan) with 
individual broad line (red) and intermediate line (blue) components and continuum (dotted 
line) shown. Note the expanded scale used to display the 2002 spectrum (middle panel). All 
fitted emission lines are listed in Table [21 
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Fig. 4. — Lya emission line profile in 2002 and 2003 STIS spectra (cyan) with individual 
broad line (red) and intermediate line (blue) components and continuum (dotted line) shown. 
The narrow velocity component of the fit (FWHM=697 km s~^, fixed to the 2002 value) is 
not shown separately here. 
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Fig. 5. — N V emission line profile in 2002 and 2003 STIS spectra (cyan) with individual 
broad line (red) and intermediate line (blue) components and continuum (dotted line) shown. 
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Fig. 6. — C IV emission line profile in 2002 and 2003 STIS spectra (cyan) with individual 
broad line (red) and intermediate line (blue) components and continuum (dotted line) shown. 
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Fig. 7. — Global CNO doublet fit from Gabel et al. (2005) to line and continuum covering 
fractions (top panel) and O VI, N V, and C IV column densities (bottom panel). Positions 
of Components 2 and 4 are labeled at top of figure. 
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Fig. 8. — Normalized O VI doublet over all observation epochs for Cjlr = Cblr = C*cont = 1- 
Positions of Components 2 and 4 are labeled at top. Top panel shows the blue doublet 
components in the 1999, 2002, and 2003 epochs together, while the second panel shows this 
for the red components. Bottom three panels show the blue and red doublet components 
together in each of the three individual observation epochs. 
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Fig. 9. — Normalized N V doublet over all observation epochs for Cjlr = Cblr = C'cont = 1- 
Positions of Components 2 and 4 are labeled at top. Top panel shows the blue doublet 
components in both the 2002 and 2003 epochs together, while the second panel shows this 
for the red components. Bottom three panels show the blue and red doublet components 
together in each individual observation epoch. 
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Fig. 10. — Normalized C IV doublet over all observation epochs for Cjlr = Cblr = C'cont = 1- 
Positions of Components 2 and 4 are labeled at top. Top panel shows the blue doublet 
components in both the 2002 and 2003 epochs together, while the second panel shows this 
for the red components. Bottom three panels show the blue and red doublet components 
together in each individual observation epoch. 
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Fig. 11. — Normalized Lyman series over all observation epochs for Cilr = Cblr = C'cont = 
1. Positions of Components 2 and 4 are labeled at top. Each panel shows an individual 
component of the Lyman series in all observation epochs: 2002 and 2003 only for Lya; 1999, 
2002, and 2003 for Ly/3 and Ly7; and 1999 and 2003 only for LyS and Lye. 
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Fig. 12.— Spectral energy distribution of Mrk 279 in 2002 and 2003, based on FUSE and 
Chandra observations. 
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Fig. 13.— Ratio of the column density of C III, C IV, N V, and O VI in the XSTAR simu- 
lations using SED2002 to those using SED2003 as a function of the 2002 ionization paramter. 
The total input column density for all is logiV = 20. 
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Fig. 14. — Blue and red components of O VI profile from 2002 FUSE data (lower left); 
blue component is repeated in other windows with profiles expected from O VI column 
density measured from 2003 FUSE data and different ILR covering fractions: Cilr = 
(red), Cilr = O.SCblr (blue), Cilr = Crlr (green). Positions of Components 2 and 4 are 
labeled at top. 
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Fig. 15. — Blue and red components of N V profile from 2002 STIS data (lower left); blue 
component is repeated in other windows with profiles expected from N V column density 
measured from 2003 STIS data and different ILR covering fractions: Cjlr = (red), Cilr = 
O.SCblr (blue), Cilr = Cblr (green). Positions of Components 2 and 4 are labeled at top. 
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Fig. 16. — Blue and red components of C IV profile from 2002 STIS data (lower left); 
blue component is repeated in other windows with profiles expected from C IV column 
density measured from 2003 STIS data and different ILR covering fractions: Cilr = (red), 
CiLR = O.SCblr (blue), Cilr = Cblr (green). Positions of Components 2 and 4 are labeled 
at top. 
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Fig. 17.— The C III A977 profile in the 1999, 2002, and 2003 FUSE observations. 
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Fig. 18. — Absorption profiles of interstellar species S III and Ar I in the 1999, 2002, and 
2003 epochs of FUSE data. 



